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A Higgs boson mass ^126 GeV as determined by the LHC data requires a large loop correction 
which in turn implies a large sfermion mass. Implication of this result for the stability of the proton 
in supersymmetric grand unified theories is examined including other experiments constraints along 
with the most recent result on cold dark matter from Planck. It is shown that over the allowed 
parameter space of supergravity unified models, proton lifetime is highly sensitive to the Higgs 
boson mass and a few GeV shift in its mass can change the proton decay lifetime for the mode 
p vK^ by as much as two orders of magnitude or more. An analysis is also given on the nature 
of radiative breaking of the electroweak symmetry in view of the high Higgs boson, and it is shown 
that most of the parameter space of universal and non-universal supergravity unified models lies 
on the Hyperbolic Branch of radiative breaking of the electroweak symmetry, while the Ellipsoidal 
Branch and the Focal Point regions are highly depleted and contain only a very small region of 
the allowed parameter space. Also discussed are the naturalness criteria when the proton stability 
constraints along with the electroweak symmetry breaking are considered together. It is shown 
that under the assumed naturalness criteria the overall fine tuning is improved for larger values of 
the scalar mass with the inclusion of the proton stability constraint. Thus the naturalness criteria 
including proton stability along with electroweak symmetry breaking constraints tend to favor the 
weak scale of SUSY in the several TeV region. Implications for the discovery of supersymmetry in 
view of the high Higgs mass are briefly discussed. 

Keywords: Higgs boson, proton decay, naturalness, LHC, Planck data 



I. INTRODUCTION 



Over the past year the ATLAS and the CMS Collaborations have identified a signal for a boson around ~ 126 GeV. 
Thus the ATLAS Collaboration finds a signal at 126.0 ± 0.4(stat) zb 0.4(sys) GeV which is at the 5.0cr level P while 
the CMS Collaboration finds a signal at 125.3 ± 0.4(stat) ± 0.5(sys) GeV at the 5.0cr level [2]. While the properties 
of the new boson still need to be fully established, it is widely believed that the discovered boson is indeed the Higgs 
boson [3]-[5] that enters in the breaking of the electroweak symmetry of the Standard Model [6l [7] . Remarkably the 
Higgs boson mass lies close to the upper limit predicted in supergravity grand unified models [SHTlj which predict 
an upper limit of around 130 GeV p!2UT6] (For a recent review of Higgs and supersymmetry see [H]). The high 
mass ~ 126 GeV requires a large loop correction which in turn implies that some of the sparticles entering the loop 
corrections (for a review see [18 ) to the Higgs mass must be in the several TeV range. In this case the heavy particles 
could be out of reach of the LHC. One possibility is that a part of the Higgs boson arises from sources outside of the 
MSSM such as from corrections arising from vector like multiplets p!9] - [22] . However, in this work we do not make 
that assumption. 

In the early analyses using radiative breaking of the electroweak symmetry (for a review see [23 ) only the Ellipsoidal 
Branch was known, in that a fixed value of the ji (the Higgs mixing parameter) implied upper limits on sparticle 
masses. However, the situation changed drastically with the discovery of the Hyperbolic Branch [24, 25 (for related 
work see [26, 27 ) when it was discovered that another branch of radiative breaking of the electroweak symmetry 
existed where the sparticle masses could lie in the several TeV region while /i could still be at the sub TeV scale. 
Specifically on this branch TeV size scalars can exist consistent with small /i. In this work we investigate the allowed 
parameter space of supergravity models under the constraint that the models accommodate the high Higgs mass. We 
show that for supergravity models most of the allowed parameter space under the high Higgs mass restriction lies 
on the Hyperbolic Branch while the Ellipsoidal Branch and Focal Point region accommodate only a small fraction 
of the allowed parameter space. We discuss the above for supergravity models with universal boundary conditions 
(mSUGRA/CMSSM) as well as supergravity models with non-universal gaugino masses (NuSUGRA). Sensitivity of 
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the proton lifetime to the Higgs boson mass is investigated and it is shown that the proton hfetime is correlated very 
sensitively to the Higgs boson mass. Further we discuss issues of naturalness in view of the large Higgs boson mass 
and the stability of the proton. It is shown that a composite fine tuning including proton stability along with the 
radiative electroweak symmetry breaking constraint prefers a SUSY scale in the several TeV region. 

The outline of the rest of the paper is as follows: In Sec.(|ll|) we discuss the radiative breaking of the electroweak 
symmetry under the constraint of the high Higgs boson mass. Here we show that most of the parameter space of 
supergravity unified models with universal boundary conditions lies on the Hyperbolic branch while the Ellipsoidal 



Branch and the Focal Point region are essentially empty. In Sec. (HI) we discuss the implications of the high Higgs 
boson mass on the proton lifetime and show that the proton lifetime is very sensitive to small shifts in the Higgs 
boson mass. Thus a shift of a few GeV of the light Higgs boson mass can change the proton lifetime by as much 



as two orders of magnitude or more. In Sec. (IV) we extend the discussion to supergravity unified models with non- 
universalities and show that the broad conclusions drawn in the previous sections still hold. In Sec.(|V|) we discuss 
the issue of naturalness and fine tuning when the proton stability constraints are combined with the constraints from 
electroweak symmetry breaking. Here it is shown that the fine tuning criteria including both the proton stability and 



the electroweak symmetry breaking constraints favor a high sfermion scale. Conclusions are given in Sec. (VI). 



II. HIGGS MASS AND BRANCHES OF RADIATIVE BREAKING OF THE ELECTROWEAK 

SYMMETRY 

It is of interest to investigate the allowed parameter space of the supergravity unified models under the constraint of 
the high Higgs boson mass. We consider first supergravity unified models with universal boundary conditions consisting 
of the universal scalar mass mo, universal gaugino mass mi/2, universal trilinear coupling Aq, tan /3 =< H2 > / < Hi > 
where H2 gives mass to the up quarks and Hi gives mass to the down quarks and leptons, and the Higgs mixing 
parameter ja which enters the superpotential via the term 11H1H2. Of specific interest is to determine the branch of 
radiative breaking of the electroweak symmetry preferred by the high mass. Thus the radiative electroweak symmetry 
breaking can be exhibited in the following form [24l [28j 

+ \mI = mlCi + A^C2 + mlj^C'^ + A/i^p , (1) 

where A'^ = Aq^ ^"^1/2 

^3 1 2^ 7 " 

tan^ p - I ^ 

Here e./, ^, k are as defined in [29^ and I^o(^) is defined by 

I)oW = (l + 6Font))-^ . (4) 

In the above Yq = /i^(0)^/(47r^), where /it(0) is the top Yukawa coupling at the GUT scale, Mq ^2x 10^^ GeV. F{t) is 
defined by F{t) = E{t')dt' , where E{t) = (1 + p^tf^^^^' (1 + ^tf^^^ (1 + pitf^^^^\ Here ft = a^(0)6i/(47r) and 
hi = (—3, 1, 11) for SU{3)^ SU{2) and U{1) and t = ln (M^/Q^) where Q is the renormalization group point. We are 
using the normalizations where 0^3(0) = 0^2 (0) = |q^i(0) = adO) and ^^(O) is the common value of the normalized 
a^s at the GUT scale. Finally, A/i^^^^p is the loop correction [30 . To understand the origin of the branches of radiative 
breaking it is useful to choose a renormalization group scale Q where the loop correction Aja'^^^^ is minimized. In this 
circumstance if all the coefficients Ci,C2,C3 are positive, the right hand side of Eq.Q is a positive sum of squares 
which leads to an upper limit on each of soft parameters determined by the size of -\- \M'^ on the left hand side. 
This is the so called Ellipsoidal Branch (EB) where ji sets an upper limit on the soft parameters and thus on the 
size of the sparticle masses. This is typically the case if the loop correction Aji^^^^ is small. However, the situation 
changes drastically if the loop correction Aji^^^^ is large. This is so because Ci are functions of the renormalization 
group (KG) scale Q and for the case when the loop correction Aji^^^^ is large the KG dependence of Ci can become 
significant. Indeed as we change the renormalization group scale Q, there is a rapid change in A/i^^^^ , and a rapid 
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compensating change also in the remaining terms on the right hand side of Eq.Q so that /i^ does not exhibit any 
rapid dependence on Q. Now it turns out that there are regions of the parameter space where one or more of the Ci 
may turn negative as Q varies. For the supergravity unified models with universal boundary conditions this is the 
case for Ci, i.e., in certain regions of the parameter space Ci can turn negative while the remainder on the right hand 
side of Eq.Q remains positive. In this case it is useful to write Eq.Q in the following form 

/+1 (EB)\ 

/.2= (FP) m^|Cl| + A^ (5) 
V-1 (HB)/ 

where stands for the rest of the terms in Eq.Q. In Eq.([5| +1 corresponds to the Ellipsoidal Branch (EB), — 1 
corresponds to the Hyperbolic Branch (HB) and Ci = is the boundary point between the two which we call Focal 
Point (FP). Its approximate form when tan/3 >> 1 is the Focus Point [31 . Ci = is achieved when Dq = 1/3 (see 
Appendix A ). We wish now to identify the allowed regions of the mSUGRA parameter space in terms of the branch 
on which they reside, i.e., EB, HB or FP. To quantify the region FP we define a small corridor around Ci = 0. This 
is feasible since FP is very sensitive to the top quark mass and we utilize the error in the top quark mass to define 
the corridor around Ci — 0. Currently the top quark mass is determined to be mt = (173.5 ± 1.0) GeV and thus we 
define the FP corridor so that |32j, 

\Ci\<5{Q,mt), ^(Q,m,)<Cl, (6) 

where 

S{Q,mt)^3{l-Do)^-^ . (7) 
mt 

and where I^o is defined in Eq.Q. Thus the Focal Point corresponds to the corridor — 1^| < Ci < |^|, the EB 
corresponds to Ci > \S\ and HB corresponds to Ci < —\S\. EB consists of closed elliptical curves and closed surfaces 
in the soft parameter space for fixed /i, while the HB region Ci < —\5\ consists of open curves and open surfaces. We 
now define a focal curve (FC) on HB as the one where two soft parameters can get large while /i remains fixed. It was 
shown in [32 that in mSUGRA there exist two varieties of Focal Curves FCl and FC2 as shown in Table I. On FCl, 
mi/2 and /i remain fixed while mo and mi/2 get large, and thus FCl is an open curve lying in the mo — Aq plane. 
On FC2, Aq and /i remain fixed while mo and Aq get large, and thus FC2 is an open curve lying in the mo — mi/2 
plane. A convolution of focal curves leads to focal surfaces ^. It is interesting to classify the allowed parameter space 
of mSUGRA in terms of the branch of radiative breaking of the electroweak symmetry they lie on, i.e., EB, HB or 
FP. This is done under the constraints of the most recent LHC searches [33^-^7 and other experimental constraints 
including the most recent results from the Planck experiment [38 . 



Focal Curve 


Large soft parameters 


Small soft parameters 


HB/FCl 
HB/FC2 


mo - Ao 
mo - mi/2 


mi/2 
Ao 



TABLE I: Classification of focal curves in mSUGRA. The focal curve HB/FCl corresponds to the case when mi/2 is kept 
fixed while mo and Ao get large keeping /j fixed (The asymptotic form of these curves give mo/Ao = ±1 27 ). The focal curve 
HB/FC2 corresponds to the case when Ao and /i are kept fixed while mo and mi/2 get large. 



We investigate the issue of classification of the branches of mSUGRA by mapping the soft parameters space in 
the following ranges: mo G (200 GeV, 30 TeV), mi/2 G (100 GeV, 5 TeV), Aq e (-6mo,6mo), and tan/3 G (1,60). 
Experimental constraints are then applied for all model points including the limits on sparticle masses from LEP [39] : 
mf, > 81.9 GeV, m~± > 103.5 GeV, m^^ > 95.7 GeV, m^^ > 89 GeV, mg^ > 107 GeV, m^^ > 94 GeV, 
and mg > 308 GeV. The most recent Planck measurement [38 of the relic density of cold dark matter gives 
ft^h'^ = 0.1199 ± 0.0027. Here we apply the 4cr upper bound, i.e. Q^h'^ < 0.13. Other constraints applied include 



The classification of the parameter space of SUGRA models into focal curves and focal surfaces is a geometric one independent of issues 
of fine tuning. The focal curves and focal surfaces automatically arise on HB for the mSUGRA case when Ci < 0. For NuSUGRA, 
the HB gets redefined such that /i remains constant while two or more soft parameters get large due to one or more of the Ci turning 
negative as discussed in Sec. IV. 
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FIG. 1: The parameter points in the mo — mi/ 2 plane in supergravity unified models with universal boundary conditions 
passing the general constraints. The plot exhibits the parameter points that lie on HB (green), on EB (red) and on FP (blue). 
The analysis shows that most of the allowed parameter space lies on HB while the allowed regions of EB and FP are essentially 
empty except for a few scattered points. For the analysis here and elsewhere in the paper we have used a top mass of 172.9 
GeV. The region excluded by the ATLAS and CMS collaborations is also exhibited. 



the constraint (-11.4 x 10"^^) < 5{g^-2) < (9.4 x 10~^) and the FCNC constraint from B-physics mea- 

surements iMS2!, i.e. (2.77 xlO-"^) < Sr (6 ^ 57) < (4.37 x IQ-"^) and Sr (5^ ^ /i+/i-) < 1.1 x 10"^ As 
done in [43] , we will refer to these constraints as the general constraints. These constraints are imposed using 
micrOMEGAs [44 for the relic density as well as for the indirect constraints and SoftSUSY [45 for the sparticle 
mass spectrum. We will also consider NuSUGRA models (for recent works on NuSUGRA see [46r48j and for 
a review see [49 . String based models also allow for non-universalities of gaugino masses, see, e.g., [5^). The 
supergravity grand unification formalism of [8 still applies. For the NuSUGRA case to be discussed in Sec. IV all of 
the experimental constraints discussed above still apply except that the ranges of the soft parameters are chosen as 



follows: mo G (200 GeV, 30 TeV), rui G (100 GeV, 5 TeV), Aq e (-6mo,6mo), tan/3 G (1,60) where z = 1,2,3 for 



In FigjTJwe exhibit the allowed parameter space of the supergravity unified models with universal boundary condi- 
tions in the mo — mi/2 plane consistent with all the constraints discussed above. The region excluded by the most 
recent ATLAS and CMS searches is exhibited. In this figure we also show the regions of the parameter space that lie 
on the HB, EB, and FP branches of radiative breaking of the electroweak symmetry. The figure shows that essentially 
all the parameter space of the universal supergravity unified model lies in the HB region (indicated by green points) 
and the EB region (indicated by red points) and the FP region (indicated by blue points) are essentially all empty 
except for a few scattered points (see also [32 ). 



In supersymmetric GUTs proton decay from dimension five operators depends very sensitively on the sparticle 
spectrum since the sparticle spectrum enters in the dressing loop diagrams which involve the exchange of squarks 
and sleptons, gluinos, charginos, and neutralinos[51 -55 (for recent reviews see [56 -58 ). Thus low values of sfermion 
masses can lead to too rapid a proton decay for the mode p vK^ in conflict with the current experimental limit [58] , 
i.e. 



NuSUGRA. 



III. PROTON STABILITY 



r 



.exp 



{p vK+) > 4 X lO^^yr. 



(8) 



Since a heavy Higgs boson mass in the vicinity of ~ 126 GeV implies relatively large values of sfermion masses it 
is pertinent to investigate proton stability within the constraint of the experimentally observed large Higgs boson 
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FIG. 2: An exhibition of the sensitive dependence of the proton hfetime for the decay mode p vK^ as a function of the 
Higgs boson mass for the supergravity unified model with universal boundary conditions. Parameters for curves 1-3 are as 
follows: Curve 1: mi/2 = 4207 GeV, = 20823 GeV, tan /3 = 7.3 while mo varies, /Mq — 50 here and for other curves; 

Curve 2: mi/2 = 2035 GeV, Ao = 16336 GeV, tan /3 = 8 while mo and Aq vary and Curve 3: mi/2 = 3048 GeV,Ao/mo = 
— 0.5, tan /3 = 6.5 while mo and Aq vary. 



mass. We will limit ourselves to generic SU{b) type models. Further, while chargino x^, gluino g and neutralino 
exchange diagrams all contribute to the decay width, the dominant contribution comes from the chargino exchange 
diagram and we will limit ourselves to considerations for decay with this exchange. Thus here the decay width is 
given by [59 , 

T{p^V,K^) = {^f\A\^m'C, (9) 

where M^/g is the Higgsino triplet mass and is the matrix element between the proton and the vacuum state of 
the 3 quark operator so that PpU]^ = eabc^a(3 < ^l^aL^bL^cLlP ^ where U2 is the proton spinor. The most reliable 
evaluation of I3p comes from lattice gauge calculations and is given [60] as /3p = 0.0118 GeV^. Other factors that 
appear in Eq.(|9]) have the following meaning: A contains the quark mass and CKM factors, Bi are the functions that 
describe the dressing loop diagrams, and C contains chiral Lagrangian factors which convert the Lagrangian involving 
quark fields to the effective Lagrangian involving mesons and baryons. Individually these functions are given by 

v2 



^msm,vlV2iALAs, (10) 



where ms{mc) are the strange (charm) quark mass, Vij are the CKM factors, and Al and As are the long distance 
and the short distance renormalization group suppression factors as one evolves the operators from the GUT scale 
down to the electro- weak scale and then from the electroweak scale down to IGeV [52| [611464] . and Bi are given by 

sm2(3msV^^ mcV2iV22 

where mf is the down quark mass for flavor i and rrit is the top quark mass. Here the first term in the bracket is the 
contribution from the second generation and the second term is the contribution from the third generation and P2, P3 
with values (±1) are the relative parities of the second and the third generation contributions. The functions Bji are 
the loop integrals defined by Bji = F(ui^dj^x^) + {dj e^), where 



F{ui,dj,x ) = [E cos J- sin j^f{ui,dj,Xi ) + C0S7+ sin7_/(?ii, d^-, X2 )] 

lfem^sin2J^i . . 7 J ~±\ 

-[£;sm7_ sm7+/(iiii,d^-,Xi ) -cos7_ cosj^f{uii,dj,X2 ) 



2 V2Mwsm(3 



-{Uil Ui 



(12) 
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and where / appearing in Eq.(12) is given by 



' Suif{uii,dj,Xk) 



' Suif{Ui2,dj,Xk) 



(13) 



Here the tilde quantities in the arguments are the sparticle masses, i.e., Ui are the up squark masses for flavor i and 
dj are the down squark masses for flavor j and the function / is defined by 



/(a, 6, c) 



ml — ml ml 



mf 



In(^) - {ma mc)]. 



m 



(14) 



Further in Eq.([l2| 7± = /3+ ± /3_ where sm2(3± = {ij.± m2)l\^v\ + (m ± r^^2)^]^/^ V^^^i = MH.(sin/3 ± cos/3) and 
sin2^^3 -2(A7 + /I cot (3)mt/{m'~^ - rn~^), E = I when sin 2/3 > /im2/M^ and E = -I when sin 2^ < /im2/M^. 
Finally C is given by 



-^[(l+ "^"^f^+^^ )(l-gf)]^ 



32^/2 ' 



(15) 



AT 



where are the stop masses and f-j^^D^F^.. etc are the chiral Lagrangian factors and we use the numerical 
values U = 0.131 GeV, D = 0.8, F = 0.47, mAr=0.94 GeV, mK=0.495 GeV, mB=1.15 GeV and we choose P2 = 1 
and P3 = —1. The partial decay lifetime of the proton into p VK^ mode is given by r(p ^ VK^) = h/T{p ^ VK^). 



Typically supersymmetric models give too rapid a proton decay for the mode p vK^ from dimension five 
operators [65^. One possible way out is the cancellation mechanism for the reduction of proton decay arising from 
different Higgs triplet representations at the GUT scale [66 . This is equivalent to raising the value of the effective 
Higgs triplet mass. ^7|. Specification of the GUT physics allows one to determine the effective Higgs triplet 
mass (see, e.g., [66^, 168]). Here, however, we do not commit to a specific GUT structure but rather consider SU{b) 
like models where due to various Higgs representations that enter at the GUT scale one has a number of Higgs 
triplets/anti-triplets Hi^Hi. Suppose we choose the basis in which only Hi, Hi couple to matter, i.e., one has 
couplings of the type [67 HiJ-\- KHi + HiMijHj, where J and K are bilinear in matter fields and Mij is the super- 
heavy Higgs mass matrix. Many grand unified models automatically lead to such a possibility [69l [70] . Specifically 
in models of the type discussed in [69^ one has only one light doublet and several Higgs triplets/anti-triplets. On 
eliminating the superheavy fields one finds that the effective proton decay operator is of the form —K{M'jj^)~^ J 
where = (M{^^)~^. This allows to be much larger than the GUT scale. In the analysis here we will 

use the effective mass M^J = M'^^/ AlAs and we consider three cases M^^ /Mq = 10, 25, 50 for analysis in this work. 

In Fig.([2| we exhibit the dependence of the proton lifetime for the decay mode p VK^ as a function of the Higgs 
boson mass under the constraints discussed in the caption of Fig.([2|. The curves show a very sharp dependence of the 
proton lifetime on the Higgs boson mass which increases by up to two orders of magnitude with a shift in the mass 
of the Higgs boson in the range of 5-10 GeV. In Fig.([3| we exhibit the proton lifetime for the decay mode p VK^ 
as a function of mo for the three values of M^J when all the parameters in the model are allowed to vary consistent 
with the radiative electroweak symmetry breaking constraints and the experimental constraints including those from 
the LHC and the Planck experiment. One finds that the parameters compatible with all the constraints clearly prefer 
values of mo in the several TeV region. 



IV. NUSUGRA: FOCAL CURVES AND FOCAL SURFACES 



In Sec. (II) a classification of radiative breaking of the electroweak symmetry is given in terms of the branches on 
which the allowed parameter space of mSUGRA resides. Here we extend the analysis to NuSUGRA and classify 
the allowed parameter space under the constraints of radiative breaking of the electroweak symmetry and all the 
experimental constraints, i.e., we discuss the composition of the parameter space in terms of HB, EP and FP. We will 
also discuss the sensitivity of the proton decay lifetime for the mode p vK^ for the NuSUGRA case. For specificity 
we define the gaugino masses at the grand unification scale by m^ where m^ = mi/2(l + ^^i), i=l72,3 and where 5i 
define the non-universalities in the U {!) , SU {2) l , SU {'iVj sectors. It is shown in Appendix [B| that in this case the 
radiative electroweak symmetry breaking equation Eq.Q for the universal soft breaking case is replaced by 

11^ + = Ciml + C2AI + C'Jmimj + Qm^Ao + A/i^ (16) 



""^ 5 10 15 20 25 30 5 10 15 20 25 30 5 10 15 20 25 30 

JH|, (TeV) mn (TeV) (TeV) 

FIG. 3: An exhibition of the partial Ufetime for the decay mode p VK^ given by blue squares as a function of mo over the 
parameter space of the supergravity model with universal boundary conditions over the allowed ranges consistent with all the 
experimental constraints. Left panel: The case when /Mg — 10. Middle panel: Same as the left panel except for the case 

M^^//Mg = 25. Right panel: Same as the left panel except for the case M^J /Mg = 50. The current experimental lower limit 
for this mode is given by the horizontal black line. The analysis given here is consistent with the Higgs boson mass within a 
2cr range. 



where Ci and C2 are as defined by Eq.([2]) while C^^ and CI are given by 



tan^/?-l ' * tan^/?- 



Here Mmn^^ and Mj are defined in Appendix B. Cs and C4, in Eq.([l7| reduce to the universal case when 
rrii = mi/2 and in this case one has C3 = -^-^ ^ 3 ^4 ~ Xli=i 2 3 ^1- -'-^ -^^S- (9) we display the dependence 

of C' s on the RG scale Q. Here one finds that in addition to Ci, and assume negative values which gives 
the possibility of new focal curves. We discuss these possibilities in further detail below. 

To examine the focal curves and focal surfaces for NuSUGRA, it is useful to define 

^3 ^1/2 — Cl^^iTTij^ = C\mi . (18) 

Further, in order to classify various regions of the radiative electroweak symmetry breaking (REWSB) for the 



NuSUGRA case it useful to write the REWSB constraint Eq.(16) in the form 



with 

3 3 

Aq (Ao + ^aimi)^, rui = ^aijmj, (20) 
i=i j=i 

where and aij are co-efiicients of linear combinations and they are functions of C2 , C^^ and CI . 

A display of the renormalization group evolution of the Ci is given in Fig. (|9| in Appendix B. Here we find that 
in addition to Ci, the elements C^^ and Cp are negative, which allows for the possibility of new focal curves and 
focal sufaces over the ones discussed in Sec.(H). Using the results of Appendix A and B one finds that four types of 
focal curves arise for the NuSUGRA case, FC1-FC4, which are listed in Table H. FCl is defined similar to the case 
for mSUGRA. FC2 has three variations: These are HB/FC2^^ where Ci > 0, C^^ < and mo and mi get large while 
Ao,m2,m3 and tan/3 remain fixed; HB/FC2^^ where Ci > 0,(71^ < and mo and m2 get large while 74o,mi,m3 
and tan/3 remain fixed, and HB/FC2^^ where Ci < 0,(71^ > and mo and m3 get large while 74o,mi,m2 and tan/3 
remain fixed. It is convenient to use the parametrization of Eq.([l8| to exhibit the effect of non-universality on focal 
curves FC2. Thus here one finds that the asymptotic value of mi/2/^0 for fixed /i as Aq gets large is affected by 
non- universality, i.e., one gets 

(21) 
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Focal Curve 


large soft parameters 


small soft parameters 


HB/FCl 
HB/FC2°^ 
HB/FC2°2 
HB/FC2°^ 
HB/FC3^^ 
HB/FC3^^ 
HB/FC4^ 
HB/FC42 


mo - Ao 
mo — mi 

mo — 1712 

mo - ms 
mi — ms 
m2 — ms 
^0 - mi 
^0 - m2 


mi,m2,ms 
Ao,m2,ms 
Ao,mi,ms 

Ao, mi, 7712 

mo, Ao,m2 
mo, Ao,mi 
mo,m2,ms 
mo,mi,ms 



TABLE II: Classification of focal curves in NuSUGRA models. Here one has the possibility of several focal curves. The focal 
curve HB/FCl is defined similar to the mSUGRA case except that mi,m2,ms are all kept fixed. As in mSUGRA here too 
mo and Aq can get large while /i remains fixed. The focal curve HB/FC2 splits into three sub cases because of the gaugino 
non-universalities. Thus the case HB/FC2°^ corresponds to the case when Ao,m2,ms are kept fixed while mo and mi can get 
large. The focal curves HB/FC2°^ and HB/FC°^ are similarly defined. For the NuSUGRA case 4 new type of focal curves 
arise. These are HB/FC3^^ HB/FC3^^ HB/FC4\ HB/FC4^ Their definitions are obvious from the table. 



>3 
1 



6^=0,6^=0 and 6^=0 
6^=0.5, Sj=0 and 5^=0 
6^ =-0.5, 5^=0 and 5^=0 
S,=0. Sj^O-S and 5j=0 
6^=0, 6^=-0.5 and 8^=0 
S.,=0, 8^=0 and 5^=0.25 
6 =0, 6 =0 and 5 =-0.2! 



/i = (0.465 ± 0.035) TeV 




5i = 
h -(-1,1) 

^3 = 

mi/2 = 2 TeV 
^0 = 1.5 TeV 

tan 13 = 45 



4 6 

mo (TeV) 



= ( 0.465 ± 0.035 ) TeV 



#{2.39,2.32 


7(1.02.0.89) 






/ HB/FC2 '"''^-^ 






I /^0.37, 0.441; 






^^^^ 








*l =-0,9,^2 = 0^*3=0 
il= 0.6,^2 = 0,^3 = 
. «i =0,42 =-0.9, *3 =0 
. «i= 0,42 = 0.9,153 = 

:|=S:l=S;|=S.i 




"^^✓^-"""''^ tan /3 = 45 

= (0.465 ± 0.035) TeV 




mo (TeV) 



FIG. 4: Top left panel: Exhibition of the Focal Curve HB/FCl of Table II with non-universalities in the gaugino sector. Here 
and in the right panel tan/3 = 45 with \i — (0.465 ± 0.035) TeV. The plot shows that non-universalities in the gaugino sector 
do not affect the asymptotic behavior of Ao/mo which is unchanged from the mSUGRA case. Top right panel: Exhibition 
of the effect of non- universalities on focal curves FC2. The analysis shows that the non-universalites have a very significant 
effect of FC2 type focal curves. The asymptotic form of the FC2 curves with non-universalities fits well with the result of 
Eq.(21). Bottom panels show the three variety of FC2 curves; left panel: An exhibition of the Focal Curve HB/FC2°^ in the 
mo — ms plane when mi = ms = mi/2 = 2 TeV and Ao = 1.5 TeV; middle panel: A display of the Focal Curve HB/FC2°^ in 
the 1710— 1712 plane when mi = ms = mi/2 = 2 TeV and Ao = 1.5 TeV; right panel: An exhibition of the Focal Curve HB/FC3^^ 
in the m2 — ms plane when mi = mi/2 = 2 TeV, mo = 1 TeV and |Ao/mo| < 0.1. The model points are colored by /i value in 
units of TeV. 



An illustration of the dependence of 7771/2/^0 on non-universalities for FC2 will be exhibited shortly. 

The focal curves FC3 arise when two of the gaugino masses get large while other soft parameters remain fixed. 
There are two possibilities here. The first one is HB/FCS''^^ where mi and 7773 get large while Ao,mo,m2 and tan/3 
remain fixed. This can happen when Ci > but C^^ is negative. The second possibility is HB/FC3^^ where 1712 and 
7713 get large while Ao,mo,mi and tan/3 remain fixed. This can happen when Ci > but Cf^ is negative. The focal 
curves FC4 arise when Aq and one of the gaugino masses get large while the remaining soft parameters remain fixed. 
There are two possibilities here. The first one is HB/FC4^ where Aq and mi get large while 7770,^2,^3 and tan/3 
remain fixed. This can happen since C2 > but C^^ is negative. The second possibility is HB/FC4^ where Aq and 
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FIG. 5: Exhibition of HB (red), EB (blue), FP (green) parameter points for NuSUGRA using the inputs given in Sec.(|n|. All 
parameter points satisfy the general constraints along with a 2cr constraint on the Higgs boson mass. As in the supergravity 
unified models with universal boundary conditions here too one finds that most of the allowed parameter space lies on the HB 
branch while EB and FP regions are highly depleted. 



1712 get large while mo, mi, ma and tan/3 remain fixed. This can happen when C2 > but Cf^ is negative. We note 
that HB/FC3^^ does not materialize since C^^ and Cf^ are both negative. Similarly HB/FC4^ does not occur since 
C2 and Cf^ are both positive. Further, while in principle HB/FC2^^ can occur when Ci < and CP is positive, 
the numerical sizes do not favor appearance of this branch. Thus as shown in the figures in Appendix pi Cf^ satisfy 
\Cl^\ <C \Cp\ <C ICf^l, where each step is roughly a factor of 10. Thus in practice the focal curve HB/FC2°^ does 
not materialize. Further, for any value of tan /3, the coefficient Ci begins positive and for tan ^5 < 5 it never be- 
comes negative (for Q < 10 TeV). Because of the above additional possibilities such as HB/FC3^^ etc are not realized. 




FIG. 6: An exhibition of the dependence of the proton lifetime for the decay mode p vK^ as a function of the Higgs 
boson mass for NuSUGRA. Parameters for the curves labelled 1-3 in the legend are as follows: Curve 1: mi = 4230 GeV, 
7712 = 843 GeV, = 3285 GeV, Aq = -27545 GeV, tan/3 = 5.3 while mo varies and M]l^ jUc = 50 here and for other 
curves; Curve 2: mi = 4794 GeV, = 3837 GeV, ms = 3856 GeV, Ao/mo = 0.842, tan/3 = 7.0 while mo and vary; and 
Curve 3: mi = 3894 GeV, = 1056 GeV, ms = 2345 GeV, Ao/mo = 2.199, tan/3 = 55.2 while mo and A^ vary. As for 
the case of the supergravity unified models with universal boundary conditions, here too one finds that the proton lifetime is 
a very sensitive function of the Higgs boson mass. 



For NuSUGRA we give a numerical illustration of some of the focal curves in Fig. Q. The left panel of the top 
row in Fig. Q gives an analysis of the Focal Curve FCl in the mo — plane. Here one finds that mo and can 
get as large as 10 TeV while \i lies in the range (0.465 ± 0.035) TeV when tan/3 = 45 and mi/2 = 0.5 TeV. We note 
that the ratio ^do/mo asymptotes to the same value irrespective of the non-universalities. A similar analysis for FC2 
is given in the right panel of Fig. Q in the mo — mi/2 plane for tan/3 = 45. Again a variety of non-universalities 
are discussed. One finds that while mo and mi/2 can get very large, i.e., as large as 10 TeV for mo and 5 TeV for 
mi/2, one still has a small /i, i.e., a \i range (0.465 ± 0.035) TeV. An analysis for FC3 is given in the three panels 
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mo (TeV) mo (TeV) nio (TeV) 



FIG. 7: Left panel: The proton lifetime shown as blue squares over the allowed parameter space in NuSUGRA models which 
pass the general constraints along with a 2a constraint for the Higgs boson mass as discussed in the text when /Mq — 10. 

Middle panel: Same as the left panel except that M^^ /Mq — 25. Right panel: Same as the left panel except that M^^ /Mq — 
50. The current experimental lower limit on the proton lifetime is shown as a black horizontal line. The analysis given here is 
consistent with the Higgs boson mass within a 2cr range. 



of the bottom row in Fig. (|4|. The left panel gives a display of the focal curve FC3^^ in the mo — m2 plane for the 
case when tan/3 = 45, Aq = 1.5 TeV, mi/2 = 2TeV and Si = = Ss and S2 lies in the range (—1,1). One finds 
that jii lies in the narrow range (0.465 ± 0.035) TeV. A very similar analysis in the mo — ma plane is given in the 
middle panel in Fig. Q where 5i = = 82 and ^3 lies in the range (—1, 1) while all other parameters are as in the 
left panel. This is the focal curve FC3^^. Finally the right panel gives an analysis of the focal curve FC3^^ in the 
Tn2 — ms plane for the case when mo = 1 TeV, mi/2 = 2 TeV, tan^^ = 45 and 60 = 0, ^2 = ( — 1, 1), and ^3 = (—1, 1). 
Here again one finds that /i lies in the range (0.465 ± 0.0350) TeV while m2,m3 get large. From a convolution 
of the focal curves one can generate focal surfaces where more than two soft parameters can vary while /i remains fixed. 

In Fig.([5| we display the nature of radiative breaking of the electroweak symmetry for all the model points within 
the allowed ranges of the parameter space for NuSUGRA. The points in red are those that lie on HB, the points 
in blue lie on EB and the points in green lie in the FP region as defined by Eqs.(|6| and As in the mSUGRA 
case here too one finds that most of the parameter points lie on HB and only a small fraction lie on EB and FP. 
In Fig. (6) we give an analysis of the sensitivity of the proton lifetime to the Higgs boson mass for NuSUGRA. As 
in the mSUGRA case here too one finds that the proton lifetime is very sensitive to the Higgs boson mass with the 
proton lifetime changing by over two orders of magnitude with a shift in the Higgs boson mass in the range of 5-10 
GeV. In Fig.Q an analysis of the proton lifetime for the mode p vK^ is given over the allowed parameter space 
of NuSUGRA within the assumed limits. The figure shows the dispersion in the proton lifetime as all the parameter 
points are varied but does show the general trend that p vK^ lifetime increases with a larger SUSY scale. 



V. NATURALNESS 



The criteria used for quantifying what is naturalness are rather subjective and various variants abound see, e.g., [161 
[241 I7TH82] . Here we discuss the fine tuning within a GUT framework including both radiative breaking of the 
electroweak symmetry and proton stability. First we discuss fine tuning for radiative breaking of the electroweak 
symmetry which is governed by the breaking condition Eq.Q. If one views as arising from the cancellation 
between pi? term and the remainder on the right hand side, it leads to a fine tuning ^24j 



An alternate criteria for fine tuning is given by the condition [71 F'^ = {a/ f{a))f' {a) where a is the sensitive parameter 
on which the function /(a) depends. Using /(a) = and the sensitive parameter as tyi^h^ finds another fine 
tuning measure 

, 2\m]r I 

F'^\^. (23) 

We will use both F and F' in the analysis for comparison. For proton decay we will use a measure of fine tuning 
defined by 
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FIG. 8: A display of the fine tuning as defined by Eqs. |22|25| vs the scalar mass mo when M^^//Mg = 50. The upper two 
panels are for mSUGRA and the middle two panels are for tJie NuSUGRA case. The left panels are when we use the fine tuning 
of Eq.(22) and the right panels are when we use the fine tuning of Eq.(23) for the electroweak sector. The red points are the 
fine tunings values for the REWSB sector, the blue points for r{p vK^)^ and the black points are the averages of the red 
and the blue points. In the bottom panel the combined fine tuning as a function of mo is given for mSUGRA (sold line) and 
for NuSUGRA (dashed line). Here we have taken the average of the left and right panels and drawn smooth curves showing 
the rapid decrease of the fine tunings as mo increases. 



This measure gives the amount of fine tuning needed in the theory parameters to enhance the lifetime so that 
the theoretical prediction is brought just above the current experimental lower limit. If we use the very crude 
approximation on the proton lifetime, i.e., r(p vK^) C • {m^± / m-M^J )~'^ and use m~ or m^± as the sensitive 
parameters, we have F' = F' 2 = 2F„c^. Thus the two ways of defining the fine tuning differ only by a small 
numerical factor. It is also useful to define a composite fine tuning by the geometric mean of the individual ones, i.e.. 



■^^m^v ' ^^^^ 

Here our view point is similar to that of [81 (for a related work see [83]). For our case n = 2 consisting of the fine 
tuning in the radiative electroweak symmetry breaking sector and the fine tuning needed to control proton decay from 
dimension five operators. An analysis of the fine tunings as a function of mo is given in Fig.(l8| where the upper panels 
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give the analysis for the case of mSUGRA and the middle panels give the analysis for NuSUGRA, and where the left 
panels give the analysis using Eq.(22) and the right panels give the analysis using Eq.(23). The red points are the 
fine tunings for radiative electroweak symmetry breaking. The blue points give the fine tuning needed in the theory 
prediction of r{p vK^) to bring the lifetime prediction just above the experimental lower limit, and the black 
points correspond to the composite fine tuning as defined by Eq.(25). One finds that typically there is a preference 
for larger values of mo for the combined fine tuning including fine tuning from the electroweak sector and the fine 
tuning needed from proton stability. This result is more explicitly exhibited in the bottom panel of Fig.([8| which 
shows fine tuning prefers regions of larger mo when the electroweak symmetry breaking and proton stability criteria 
are combined. A similar conclusion was arrived at in the work of [81] which combined the electroweak symmetry 
breaking, FCNC and CP violation criteria. 



VI. CONCLUSIONS 



The high mass of the Higgs boson discovered recently requires a large loop correction to its mass which points 
to the possibility that the overall weak scale of supersymmetry may lie in the several TeV region and could even 
be as large as tens of TeV. If the scalar masses are that large they would help resolve one of the serious problems 
of supersymmetric grand unification related to proton decay. Thus proton decay from lepton and baryon number 
violating dimension five operators often leads to proton lifetimes which fall below the current experimental limits. In 
this work we show that the proton lifetime is a very sensitive function of the Higgs boson mass in a unified theory. 
Thus a few GeV upwards shift in the Higgs boson can result in orders of magnitude suppression of the proton 
decay from baryon and lepton number violating dimension five operators and a corresponding enhancement of the 
proton lifetime. The analysis is first done for the mSUGRA model and then extended to NuSUGRA. Here we also 
analyse the allowed parameter space in terms of which branch of the radiative breaking of the electroweak symmetry 
the parameters lie, i.e., whether on the Ellipsoidal Branch, the Hyperbolic Branch or the Focal Point region. The 
analysis presented in this work shows that under the current experimental constraints including those from LEP, 
Tevatron, LHC, FCNC and the Planck data ^84^ one finds that most of the parameter points of mSUGRA and of 
NuSUGRA models lie on the Hyperbolic Branch with only a very small fraction lying on the Ellipsoidal Branch or in 
the Focal Point region. We also discuss issues of naturalness and fine tuning and show that the composite fine tuning 
including fine tuning from the electroweak sector and from the stability of the proton points to high scalar masses. 
However, some of the gauginos can be light with their masses mostly limited by their lower experimental limits. 
These include the light chargino, the lightest neutralino, and the second lightest neutralino and the gluino. These 
should be accessible with increased energy and luminosity at the next round of experiment at the LHC. Regarding 
proton decay discovery of the supersymmetric mode p VK^ is over due and this mode continues to be the most 
likely candidate to be discovered first in the next generation of proton decay experiments. 
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Appendix A: The relation of Ci(Q) to 



Here we establish the relation between the Ci{Q) and mH2' The RG evolution connects and on the third 
generation masses mu and tuq and one has 
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Here Yi = /i^/(47r^) where hi is the top Yukawa coupling and Ai is the trilinear coupling in the top sector. The 
above equations with universal boundary conditions at the GUT scale allow a homogeneous solution satisfying [3T] 
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where J is an integration factor defined by 



J(t) =exp -6 / Yt{t')dt' . 
I Jo 



(A2) 



(A3) 



As Q ^ Mq^ one has J it) 1 and the universahty of the masses is recovered at the GUT scale. In [32 a connection 
was estabhshed between Ci{Q) and dmH2 which we now ihustrate. Thus Y{t) at the one loop level satisfies the 
equation 



dt 
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and one finds 



Ytit) = 



YiO)Eit) 
1 + QY{0)F{t) 



(A4) 



(A5) 



where F{t) and E{t) are defined after Eq. Q. It is then easy to see that J{t) = D^ii)^ where 0^(1) is defined by 
Eq. Q. Thus 5m\ takes the form 



Sttl 1 



(A6) 



and Ci can be expressed in terms of Sm'jj^ 



tan^ P 



SrUff tan' 



(A7) 



Ci = was defined as the Focal Point in [32 . At the Focal Point ji^ essentially becomes independent of mo- For 
tan/3 >> 1, Ci —^^1/2 vanishing of Ci implies vanishing of Sfnjj^ which is defined to be the Focus Point. 

Thus the Focal Point defined by Ci = is just the boundary point between EB defined by Ci > and HB define 
by Ci < 0. For the NuSUGRA models all solution where some of the soft parameters can get large while /i^ remains 
fixed lie on HB. This can happen when some of the Ci other than Ci turn negative, as discussed in the Appendix 
below. 



Appendix B: Analysis of C^s for Models with Non-universalities in the Gaugino Sector 

The presence of non-universalities in the gaugino sector affects the co-efficients Ci and in this Appendix we give a 
computation for these by inclusion of non-universalities in the SU{3)c^ SU{2)l and U{1) gaugino sectors. We begin 
with the radiative electroweak symmetry breaking with the inclusion of non-universalities in the gaugino sector. We 
have 

^ (tan/32-l) 2^z^^^. (Bl) 



with 



= ml+ (^^/i + , (B2) 
= e{t) + ^ofit)+mlh{t)-Alk{t), (B3) 



and fi{t) is defined by fi{t) = z( mf where 
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It is useful to introduce a column vector = (mi, m2, ms) and a matrix Mmn^ such that m|^^ = rfi^ • Mmn^ ' ^ = 
(^MmH-^)^j ^i^j where Mh^ is given by 
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2 / 

V 0/ 



(B5) 



Thus we have 



m 



rujirij. 



(B6) 



The above exhibits the gaugino mass dependence of m|^^ explicitly. Now let us look at m^^^ given by Eq.(B3) and 
write it in a form which exhibits the gaugino mass dependence explicitly. Now m|^^ contains the functions e{t) and 
/(t) which are given as 
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where Hi(t) are defined by 



H4 = F{t)H2{t)-H3{t), H5 
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and hi are defined by hi = Z^rrii with 



H2{t) then takes the form 



where Mir is a row vector 
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Similarly, we may write all the M^. (t) in matrix or vector forms so that Hs = • m, H4 = M^^ • m, 
mJ ' • m, H6 = (^^h2 ' ^) • ^8 = • ' ^ where the matrices etc are given by 

M^3= f E{t')M^^{t')dt', M^^=M^^{t) f E{t')dt'- f M^^{t')E{t')dt', 
Jo Jo Jo 
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FIG. 9: The upper panels: RG evolution of Ci{Q) and C2{Q) as a function of the renormalization group scale Q at different 
tan/3. Left panel: Ci(Q) at tan/3 = 5, 6, 10 and 45. Right panel: C2(Q) at tan/3 = 5, 6, 10 and 45. It is seen that Ci(Q) turns 
negative as the scale Q increases while C2{Q) remains positive. It is also seen that Ci{Q) is very sensitive to tan/3 while C2 
is very insensitive to tan/3. The lower panels: An exhibition of Cs at different tan/3. Left panel: C^'^ at tan /3 = 5,6, 10 and 
45. Middle panel: at tan/3 = 5,6, 10 and 45. Right panel: Cf^ at tan/3 = 5,6, 10 and 45. It is seen that €3^ and are 
negative, which allows the possibility of new focal curves as discussed in the text. 
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Similarly, Fi{t) defined by 



F2 = ^/1 + -/2, 
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^3 = F{t)F2{t) - I E{t')F2{t')dt' , 
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can also be written in matrix forms so that F2 = [Mp^^. rriimj^Fs = {^Mp^^. mirrij^F/^ = [Mp^^. mirrij, with Mp, 
defined by 
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-z^j 



Mp^it) = F{t)Mp^{t) - E{t')Mp^{t')dt', Mp^it) = / E{t')Mfj^{t')dt'. 
Jo Jo 

We repeat the same procedure for functions Gi defined by 

Gi = F2{t) - \H2{t)\ G2 = QF3{t) - F^{t) - 4H2{t)H4{t) + 2F{t)H2{tf - 2He{t), 



(B19) 
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(B21) 
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which could also been written as, Gi{t) = (M^J.. niimj, G2{t) = (M^,^) .. minij, with Mq, defined by 
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(B22) 
(B23) 



We return now to e{t) and f(t) and write these in the matrix form so that e(i) = {Mq)^. mirrij, and f{t) = (^Mj^ rui 
with 



1 



Using the above we can write m'jj^ in the form 

mjj^ = (Me) mirrij + Aq (^Mj^ rrii + mlh{t) - Alk{t). 



(B24) 



(B25) 



Thus using Eq.(B6) and Eq.(B25) in Eq.(Bl), we finahy have the radiative electroweak symmetry breaking equation 
for non-universaUties as given in Eq.(n!6|). 
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